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ABSTRACT 


The Louis et al. (1982) bulk aerodynamic method for air-sea flux estimates is 
currently used in mesoscale models such as COAMPS, while the TOGA-COARE method 
is a state-of-the-art flux parameterization involving recent findings in surface layer 
meteorology, and has been proposed as a replacement to Louis. The Louis and TOGA- 
COARE methods were compared to direct observations during ACE-1. 

Results firom both methods compared well to observations for momentum flux. For 
sensible and latent heat flux, both methods showed good agreement in low flux regimes 
and underestimated the fluxes at higher values. Calculation of the transfer coefficients 
required to match observations indicated that the bulk transfer coefficients do not increase 
rapidly enough for z/L values less than -0.5. In the high wind regime, the transfer 
coefficients were very sensitive to static stability. The CO ARE method was superior to the 
Louis method for sensible heat flux estimates while the Louis method was better for latent 
heat flux. A sensitivity test was done to use the COARE roughness length 
parameterization in the Louis method. This resulted in slight improvement in sensible heat 
flux estimates for highly convective conditions. Latent heat flux was overestimated by the 
modified Louis parameterization in the same manner as the COARE algori thm ^ indicating 
that specification of the latent heat roughness length requires further study. 
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I. INTRODUCTION 


Surface fluxes of momentum, heat, and water vapor 
represent the exchange of energy and mass at the interface 
between the ocean and the atmosphere. These fluxes are the 
driving force for atmospheric and oceanic boundary layer 
growth, mesoscale ocean current and wave evolution, and 
tropical cyclone genesis and intensification. As a result, 
an accurate representation of these air-sea fl\ixes is 
critical to numerical prediction of a variety of weather, 
climate, and ocean problems. The increased interest in 
coupled and high resolution air-sea numerical models has 
resulted in a pressing need for more accurate surface fliix 
parameterizations that are still computationally efficient. 
As numerical weather prediction becomes more focussed on 
mesoscale phenomena that are driven by the local forcing, 
errors in the surface flvixes result in errors in depicting 
the scales of interest and in capturing the evolution of the 
boundary layer on both sides of the interface. 

Turbulent fluxes can be calculated in a variety of ways 
depending on the availability of the measurements. The 
commonly used methods include the eddy correlation method, 
the inertial subrange dissipation method, the flux-profile 
method, and the bulk aerodynamic method. Of these, only the 
eddy correlation method calculates the turbulent flxoxes 
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directly using a time series of high frequency measurements 
of wind, temperature, and humidity (Chou, 1993). The 
advantages of this method is that it is direct and 
conceptually straight forward. Disadvantages are that 
expensive fast-response sensors must be used, the high-rate 
measurements thus obtained are not routinely available from 
weather stations or at-sea sources, and special platforms 
are required to prevent flow distortion around the sensors 
(Stull, 1988; Chou, 1993). 

The dissipation method makes use of the spectral 
characteristics of the inertial subrange. For steady state 
turbulent flow this energy cascade produces a -5/3 slope on 
a log-log plot of the spectrum, from which the dissipation 
rate of the Turbulent Kinetic Energy (TKE) can be estimated. 
The frictional velocity is then obtained from the 
dissipation rate using Monin-Obukhov similarity theory. The 
major advantage of this method is that the momentum flux 
thus obtained is not sensitive to uncertainties associated 
with moving platforms such as a ship or buoy. The method is 
thus widely used with ship and buoy measurements. The 
disadvantage is that, similar to the eddy correlation 
method, fast-response sensors must be used to produce the 
appropriate time series of perturbations (Stull, 1988). 

The flux-profile method uses the integrated form of the 
flux-profile equations based on Monin-Obukhov similarity 
theory and the empirical flux-profile relationships. In 
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general, two levels of measurements are needed, and the 
fluxes or scalar profiles are obtained through iteration. 
When only one level of measurement is available, the second 
level of information is obtained by parameterizing the 
height of the lowest altitude in the surface layer (Liu et 
al., 1979). The advantage of this method is that it can be 
applied to mean quantities from bulk measurements. The 
disadvantage is that it requires two levels of measurements, 
which are not routinely available from ships or 
meteorological stations. 

The essence of the bulk aerodynamic method is to 
parameterize surface fluxes with mean quantities routinely 
available from conventional data sources. A variety of bulk 
aerodynamic schemes have been developed by different 
researchers for various mean environmental conditions. All 
of these methods are based on Monin-Obukhov similarity 
theory which makes use of generalized properties of the 
atmospheric surface layer. The fliixes are then estimated 
based on averaged measurements at a single level. 
Theoretically these are ensemble averages but more often 
they are space or time averages. The simplest of the bulk 
methods use fixed transfer coefficients based on vertical 
integration of empirically derived flux-profile 
relationships. More sophisticated models calculate the 
transfer coefficients as a function of static stability and 
wind speed. These transfer coefficients are then multiplied 
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by the wind speed profile between the ocean current and the 
measurement level and the moisture, temperature, or speed 
profile as appropriate for the flux of interest. In addition 
to a dependence on stability, the literature has also 
addressed the dependence of heat flux transfer coefficients 
on wind speed. Wind speed and static stability greatly 
affect the accuracy of these schemes since the profiles of 
heat, momentum and moisture are dependent on the turbulence 
structure. The turbulence structure is in turn dependent on 
the buoyancy and mechanical production in the boundary layer 
and in response to the large scale environmental forcing. 

The advantage of the bulk method is that it uses routinely 
available measurements of air temperature, sea surface 
temperature, and wind speed to estimate the fluxes and is 
thus the most suitable for operational numerical weather 
prediction. The disadvantage is that significant error can 
be introduced in the model in regions of complex or 
heterogeneous conditions or when the synoptic conditions 
vary largely from the conditions in which the 
parameterization was developed. Furthermore, the bulk 
aerodynamic method is by definition an estimate of the 
fluxes based on mean conditions which may not capture fine- 
scale structure significant in mesoscale models. 

In numerical forecast models the consideration of 
accuracy with varying turbulence structure must be weighed 
against computational expense. In the past, the Louis 
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surface flux parameterization (Louis, 1979, Louis et al., 
1982) has been widely used because it includes an assessment 
of stability dependence but is a non-iterative approach 
using an estimated bulk Richardson number to achieve 
closure. This results in a compact parameterization which is 
reasonably accurate and computationally efficient. 

More recent work has established that the calculated 
transfer coefficients are too large in the stable case 
(Holtslag and Beljaars, 1989) and several modifications 
applying better empirical formulas for stable stratification 
to and to account for different roughness lengths for heat 
and momentum have recently been proposed (Beljaars and 
Holtslag, 1991; Launianen, 1995; Uno et al., 1995; and Lo, 
1996). The current parameterization used in the United 
States Navy's operational mesoscale model, the Coupled Ocean 
Atmosphere Mesoscale Prediction System (COAMPS) is based on 
the Louis et al. (1982) method, as are several other 
mesoscale forecast models such as the Regional Atmospheric 
Modeling System (RAMS) (Pielke et al., 1992) and General 
Circulation models such as the BMRC AGCM (McAvaney and Hess, 
1996). There have also been some studies indicating that 
Louis-type flux parameterizations work reasonably well in 
neutral and stable conditions but may break down in 
convective conditions (Garratt et al. 1996). 

Another popular parameterization is the Liu-Katsaros- 
Businger (LKB) scheme (Liu et al.,1979) which uses the 
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roughness Reynolds number instead of the bulk Richardson 
number to achieve closure and uses iteration to 
simultaneously solve the surface layer equations based on 
Monin-Obukhov similarity theory. See Fairall et al. (1996b), 

DeCosmo et al. (1996), and Chou (1993) for discussions of 

this method when compared to field data and proposed 
modifications to improve accuracy. 

In this study we will examine the modified bulk 
aerodynamic parameterization scheme used in the current U. 

S. Navy's mesoscale model (COAMPS) and compare the results 
to direct measurements using the eddy correlation method. 

The bulk parameterization, which is discussed in more detail 
in section II, is based on Louis (1979) but has undergone 
several modifications. The fltix parameterization developed 
during the Tropical-Ocean Global Atmosphere Coupled-Ocean 
Atmosphere Response Experiment (TOGA-COARE) (Fairall et al., 
1996b) is based on the LKB model and is specifically 
designed to produce more accurate flux estimates from bulk 
methods. This scheme will be applied to the same dataset in 
order to assess its skill over the Louis method as compared 
with direct flux measurements derived from the eddy 
correlation method. The TOGA-COARE algorithm is proposed as 
a replacement for the Louis model in COAMPS (Burk, personal 
communication) but the stability functions have been 
replaced with a polynomial fit to reduce computational 
expense. Despite the possible need for another correction to 
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roughness lengths in the wave breaking regime as discussed 
in DeCosmo et al.{1996), this model does allow for separate 
computation of the roughness lengths for heat and moisture, 
adjusts for skin temperature versus bucket temperature 
inputs to the calculation, and incorporates other 
improvements that are expected to produce more accurate 
fluxes than the original LKB scheme. 

The objective of this research is to examine the 
accuracy of the Louis surface flux parameterization as a 
widely used method in current numerical weather prediction 
as compared to the TOGA-COARE algorithms which represent the 
state-of-the-art in accuracy but are more computationally 
expensive. Direct flvix measurements calculated using the 
eddy correlation method are used as a benchmark. The dataset 
used in this study was collected during the Southern Aerosol 
Characterization Experiment (ACE-1) conducted in the Tasman 
Sea in November and December 1995 (Bates et al., 1998) and 
includes a variety of cases including stable and high wind 
situations. This database is therefore ideal for a thorough 
examination of surface flux parameterizations. 
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II. BACKGROUND 


A. THE BULK AERODYNAMIC METHOD 

In the surface layer, generally the lowest 10% of the 
boundary layer, momentum, heat, and moisture flux are nearly 
constant. These fluxes can be represented as: 

T =p/ {w'u’) {w'v') ^=puf 


H=-pCw'@'=-pcuB 

^ P p ♦ ♦ 


X=—pLyq'=-pL^u^q, 


where x is the wind stress, H is the sensible heat flux, X 


is the latent heat fliix, p is the air density, Cp is the 
isobaric specific heat, and is the latent heat of 


vaporization. The bulk aerodynamic relationships 
approximating these fluxes are given in Eqs. (l)-(3): 
T=pC {U-U )2 

jn ' s' 


( 1 ) 


H=-pCy(U-U^) 

K=-pLy^U-U^) iq^-q) ^ 

Measured quantities are the surface (microlayer) potential 
temperature ( 63 ) which is assumed to be in the 2 rmal 
equilibrium with the sea surface potential temperature ( 63 ) , 
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potential temperature (6), specific humidity (q), and 


saturation specific humidity {q^) as a function of sea 

surface temperature (SST) and sea level pressure (SLP). 

is the current velocity. It is usually small and assumed to 
be zero in moderate to high wind conditions. The transfer 
coefficients (C„ , Ct^ , Cg) can be calculated based upon 

surface layer similarity theory. Integrating the flux- 
profile relationships between the respective roughness 
lengths (zom, Zoh/Zoq) the measurement height z, we obtain 

the relationships: 

k(U-U^) 


k(e-ej 


ln(z/z„,) -¥,(z/i)+T,(z„,/L) 


( 6 ) 


where u* is the friction velocity scale, 6. is the 
temperature scale, g* is the humidity scale, and k is the 
von Karman constant (taken to be k=0.4 in Louis et al. 
(1982)) . , and Wg are the integrated forms of the 
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empirical stability functions. The subscripts m, h, and q, 

denote the quantities for momentum, sensible heat, and 
latent heat respectively. L is the Monin-Obukhov length 

defined as: 


L= 



(7) 


where g is gravity. 

The two parameterizations examined here differ in their 
approach to estimating the surface fluxes. The Louis method 
uses empirically adjusted curve fitting to approximate the 
stability functions and defines the appropriate scaling and 
stability regime by the momentum roughness length (zoJ , 

measurement height (z), and the Richardson number (R^s) . The 

TOGA-COARE algorithm calculates the roughness lengths for 
momentum, heat, and moisture based on the characteristics of 
the interfacial layer and then treats this as a second 
measurement level to solve iteratively for the surface layer 
fluxes from the flux-profile relationships. 

B. THE LOUIS PARAMETERIZATION 

Following the treatment by Louis (1979), surface fluxes 
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for momentum (i), sensible heat(H) , and latent heat (A) are 

given by Eqs. (l)-(3). The original development combined 
latent and sensible heat into a single term but these can be 
separated by retaining the same basic form of the equations 
for each parameter. 

In order to get u., 6., and g*, from measurements of the 

mean quantities at a single level, one must solve Eqs. (4)- 
(7) iteratively, particularly in unstable conditions. This 
process can be very computationally expensive, particularly 
in high resolution or large domain size models. Louis (1979) 
developed a more computationally efficient approach. He used 
the bulk Richardson number (Ri^) , given in Eq. (8) applied 

to a derived set of functions (F„, F^, and Fg) that closely 

match the iterative solution for the unstable case and field 
experiment data for the stable case. Equations (9) and (10) 
are the relationship between the fitted functions and the 
fluxes. The turbulent Prandtl number (R) in Eq. (11) is 

empirically derived from data collected over land. C„n and 
Chn are the momentum and heat transfer coefficients in 
neutral conditions. 

_ gzAe 

7r..2 



(9) 
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(10) 



The fitted relationship between fluxes, and R^s was 

incorporated into the global model at the European Centre 
for Medium-Range Weather Forecasting (ECMWF) where further 
modifications were made in order to obtain realistic model 
results. In the stable case, the fitted flux relationship 
was found to vanish at a critical Richardson number value of 
0 .21, resulting in the numerical model surface becoming 
energetically disconnected from the atmosphere. With no heat 
flux from the surface into the boundary layer, the 
atmosphere therefore cooled too fast. Based on approximating 
the analytic solution for the unstable case and field data 
in the stable case, Louis (1979) gives the explicit 
expressions of F as: 

bRi_ 

F=1 -- —, unstable 

l+clRigl^/^. 

® ( 12 ) 
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F= 


1 

n+b'Ri^)^ 


stable 


(13) 


where b, b', and c are empirically derived constants. 


Louis et al. (1982) modified this parameterization 
based on diagnostics of the global operational model at 
ECMWF and additional theoretical and consistency 
considerations. The final set of equations have the form: 


C=a^(l-- 


2bRi, 


l+2ba^c{-^\Ri^\)^'^ 


), unstable, b=5, c=l. 5 


(14) 


C. =aUl- 

n, g 


2bRi, 


), unstable,b=5, c=5 


l+3ba^c{ — \RiJ) 


1/2 


(15) 


C =a^ (-) , stable, b=5, d=5 

" 1^2bRi,{l^dRi,)^^^ 

C. =a^{ - - -—), stable, b=S,d=S 

k=0.4 (18) 

where C„,h,q=a^F„,h,q- 
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There have been some modifications to the Louis et al. 
(1982) algorithm. Meaningful observations have never been 
easy to obtain in stable conditions due to lower levels of 
turbulence and the influence of inhomogeneities in the 
terrain (Garratt et al., 1996) but carefully designed 
experiments and more recent theory have improved the 
parameterization in stable boundary layers (Beljaars and 
Holtslag, 1991). The parameterization in COAMPS has been 
modified to account for different roughness lengths for heat 
and momentum and to use updated empirical formulations for 
transfer coefficients in the stable case. The momentum 
roughness length (zom) is estimated from the u* of the 
previous time step using Charnock's relation (1955), 
different roughness lengths for heat and momentum based on a 
fixed ratio of l/lO are applied, and a modified set of 
empirical functions for the stable case have been 
substituted based on work by Holtslag and Beljaars (1989). 
The modified functions in the stable case are: 


F =F =F =■ 


l+lOHi^d+SRi^l 


•,stable 


(19) 


k=0.4. 



( 20 ) 
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C. THE T06A-C0ARE SURFACE FLUX ALGORITHM 

The TOGA-COARE surface flvix algorithm is a state-of- 
the-art method with the stated goal of flux accuracy of 10 W 
m"^. The algorithm is based on the LKB method and has been 
modified to accommodate various difficult environmental 
conditions such as rough sea and the free convection limit. 
It also includes several physical considerations such as the 
cool-skin and warm-layer correction for the SST inputs and 
roughness length effects neglected previously because their 
impact individually was considered to be small. 

Fairall et al. (1996b) provides a more detailed 

description of the algorithm which was applied to 
observations collected during TOGA-COARE. The algorithm 
calculates surface fluxes from the fliix-profile 
relationships expressed in (4)-(6), which call for the mean 
quantities at two known levels. While the observations at a 
reference height provide the inputs at one level, the values 
at a second level come from the roughness height, which is 
defined as the lowest level in the surface layer, or 
alternatively the top of the interfacial sublayer. One 
modification in the TOGA-COARE flux algorithm to the 
original LKB method is the specification of the roughness 
length . Equation (21) combines the LKB formulation for 

smooth flow (wind speed less than 2 m s'^) with Charnock's 
relation, which is valid for rough flow (wind speed greater 
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than 8 m s"^) . Appropriate heat and moisture roughness 
lengths Zpg) are more uncertain. The current TOGA- 

COARE scheme parameterizes these roughness lengths as a 
function of the roughness Reynolds number (Rr) as discussed 

in section II.D.3. Equations (23) and (24) are the neutral 
transfer coefficients at the reference height in the TOGA- 
COARE model where R is the Prandtl number (R= 1.0 in 


Fairall, 1996b). Equation (22) gives the roughness Reynolds 
number (Rr) and Eqs. (25)-(28) are the empirical stability 
functions following Fairall et al. (1996b): 

2 


V 

z._=a—+0.11-^ 


om 


( 21 ) 


R = 

r 




( 22 ) 



K 

[log(^)-4r^(a] 


(23) 



_Rk_ 

[log(^)-Hr,,^(01 

0h,q 


(24) 
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( 25 ) 


W =2*ln[ 

m ^ 


(l+X) 


] +ln [ — ^ ] -2 *arctan (X) +-^, unstable 






], unstable 


( 26 ) 


X= 




1/4 




1/2 


, ^=z/L,a=a. =16 

^ ' J7J ^ 


(27) 


’P, 


„ . =-jb„ h I, stable,b=b. =1 


(28) 


where ^=z/L where L is the Monin-Obukhov length, defined in 
Eq. (7). 

The TOGA-COARE parameterization is particularly 
concerned with accurate flux calculations for energy budget 
analysis in the equatorial Pacific. Therefore, another 
important modification to the LKB scheme is in the 
vertically integrated stability functions to include the 
free convection limit when u. approaches zero but buoyancy 

flux does not. The vertical gradients of the scalar 
quantities are expected to show a dependence which 

leads to the convective profile in Eq. (29). Equation (30) 
combines the standard Businger-type profiles in (25) and 
(26) with the free convective case in (29). 

^ =1.5 *ln [ ] -TIarctan () +—, y=12 .87 

3 ^3 ^ ( 29 ) 
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( 30 ) 


(jn,h,g) 


1 




+ —Z^ 
1 + ^2 0 


Other modifications to the LBK formulation which do not 
^pply to this study and therefore will not be discussed in 
detail include a gustiness (w^) factor which accounts for 

flux induced by boundary layer scale variability in calm 
conditions, the Webb correction (^) which accounts for the 

requirement that mean vertical motion must occur associated 
with the heat flxox so that net dry mass flux is zero, and 
the sensible heat contribution carried by precipitation 
(Hsr) • The justification for these corrections is detailed 

in Fairall et al. {1996b). These modifications affect the 
total flux estimate from the LKB formulation by about 20%. 
The TOGA-COARE algorithm has been verified in low wind 
regimes and convective conditions but may still underpredict 
the momentum fl\ix in winds above 10 m s'^ as the LKB scheme 
was found to do in De Cosmo et al. (1996). 

To further reduce the error in surface flux estimation, 
the TOGA-COARE fl\ix algorithm includes corrections to adjust 
the bulk sea surface temperature (SST) to the skin 
temperature. These include the warm layer and cool skin 
effects discussed in section II.D.2. 

The algorithm is applied as follows: 

Step 1. Input measurement height, latitude, and 
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longitude, water temperature sensor depth, inversion height, 
and surface pressure, and set all constants. 

Step 2. Input measured variables: u, 6s, T, q, and rain 

rate (J?) . Correct and for warm layer and cool skin 

effects as discussed in Section II.D.2 

Step 3. Assign W5=0.5 m s"^ and neutral transfer 

coefficients as a first guess and compute all temperature 
dependent constants and initial values for u., T., and g.. 

Step 4. Modify the functions for stability by iterating 
Eqs. (l)-(3), (21)-(24), and (28) - (30) , until convergence, 

which is normally within five iterations. 

Step 5 & 6. Compute fluxes and wind stress, increment 
integrals for the warm layer correction, and go to the main 
loop. 


D. PHYSICAL CONSIDERATIONS 

1. Wind Speed 

In the bulk parameterization scheme, the presence of 
wind speed denotes the role of the vertical wind shear in 
generating turbulent mixing and therefore turbulent fluxes. 
Since the parameterization is based on Monin-Obukhov 
similarity theory, wind conditions that violate this theory 
will result in problems with the bulk method. One of these 
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special wind regimes is in the limit of free convection when 
the mean wind speed is close to zero. Modifications were 
made in TOGA-COARE to include a gustiness wind factor in the 
apparent wind speed and a modified stability function to 
account for free convection (Fairall et al. 1996b). 

High wind conditions over the ocean also can violate 
the assumptions of similarity theory and the bulk flxix 
formulation. Based on observations during the Humidity 
Exchange Over the Sea (HEXOS) experiment, DeCosmo et al. 
(1996) indicates that in winds equal to or greater than 10- 
15 m s‘^ flux-profile relationships for latent and sensible 
heat may break down. This may be due to changes in the 
surface layer structure due to waves, wave-breaking, and 
droplet evaporation. Whitecaps begin to form at winds of 3 m 
s'^ and breaking waves cover 1% of the surface at winds 
above 10 m s"^, increasing with the ciobe of the wind speed. 
These effects may also be dependant on other variables such 
as wave age and direction in relation to the wind, which is 
not always well known. In addition to numerical modeling and 
synoptic models, flux parameterizations in high wind are 
important for better understanding and prediction of 
tropical cyclone genesis and intensification (see for 
example. Black and Holland 1995). The HEXOS observations 
reported in De Cosmo et al. (1996) indicated that the flux 
exchange coefficients do not significantly change with wind 
speed for water vapor flux up to speeds of 18 m s"^ and 
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sensible heat flux up to winds of 23 m s'^ and that the LKB 
scheme is valid at these speeds for heat and moisture. 
However, the drag coefficient and momentum roughness length 
were found to increase more rapidly than LKB predicts for 
wind speeds greater than 10 m s'\ The work by Smith (1980) 
and DeCosmo et al. (1996) indicates that bulk transfer 
coefficients for heat and moisture flux are not seen to 
change with increasing wind speed. This contradicts the 
results of Liu (1979) and Fairall et al. (1990) that show a 
decrease of the bulk transfer coefficients with increasing 
wind speed, however scatter in the data prevents a 
definitive result. As discussed in Fairall et al. (1996b) 
and DeCosmo et al. (1996), the data indicates a consensus 

error of about 30% in estimation of the bulk transfer 
coefficients so the issue can not be definitively resolved 
without more precise formulations. Since the increase in u* 

with wind is not accompanied by an increase in the heat and 
moisture exchange coefficients, the roughness lengths for 
heat and moisture must decrease with wind speed, as in the 
LKB formulation above a u. value of 0.2. 

2. Sea Surface Temperature 

Bulk parameterizations of sensible heat flux are 
particularly sensitive to accurate measurements of sea 
surface temperature (SST). Estimates of momentum flux can 
also be affected by SST accuracy since a sounding is usually 
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not available in numerical forecast models so stability is 
estimated by the air-sea temperature contrast. The 
theoretical development requires 6^ to be the skin potential 

temperature, which ideally would be the measurement of by 

an infrared radiometer (Fairall et al. 1996a). This method 
requires careful correction for reflected atmospheric 
i^sdiance and close tolerances in terms of absolute accuracy 
and drift. Inexpensive, accurate, and reliable instruments 
are not available to directly measure the surface radiance 
so in situ measurements from sensors placed in the water are 

often used. Another means of measuring ocean surface 
temperature is from satellite-derived multi-channel sea 
surface temperatures (MCSST). These values are appropriate 
for most bulk methods since they measure the average 
radiance in a area related to the satellite resolution, 
nominally 8-25 kilometers on a side for the NOAA polar 
orbiters (Kidder and Yonder Haar, 1995). Accuracy is on the 
order of 0.6 K for the operational AVHRR instrument when 
compared to ship observations. MCSST retrievals use 
regression techniques to match the raw radiances to ship and 
buoy observations. These observations should therefore be 
considered bulk temperature values and not skin 
temperatures. 

Ship and buoy based bulk sensors placed directly in the 
water generally either use floating sensors at a few 
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centimeters depth or measure the water temperature of the 
ships intake at 2-10 meter depth. As discussed in Fairall et 
al. (1996a), to obtain the correct interfacial temperature 
these measurements must be corrected for the warm layer and 
cool skin effects. The warm layer effect refers to the 
result of the fact that about half the solar radiation 
received is absorbed by the upper meter of the ocean. This 
leads to substantial diurnal variation in the profile of the 
first few meters of ocean mixed layer (Price et al. 1986). 
The temperature difference between the interfacial layer and 
the sensor inlet can be as much as several degrees in the 
first meter in very light winds with strong solar radiative 
flux. The cool skin effect refers to the fact that because 
sensible heat, latent heat, and longwave radiant fluxes 
occur in the upper fractions of a millimeter at the surface, 
a floating sensor at even a few centimeters will measure a 
warmer temperature than is at the interface. This "cool- 
skin" can be as much as 0.2 to 0.5 K lower than the water a 
millimeter below the surface (Fairall et al. 1996a). The 
ocean convective effect causes the cool skin effect to reach 
a maximum at night and with decreasing wind speeds. 

3. Roughness Length 

An implicit assumption in Louis et al. (1982) is that 

the roughness lengths for momentum (zo,,) , sensible heat {zoh) 

and latent heat (z^g) equal. This assumption has been 
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examined in field studies by Beljaars and Holtslag (1991) 
and numerical studies by van den Hurk and Holtslag (1997) 
and Lo (1996) which have shown this assumption is only valid 
in certain specialized cases and can produce significant 
errors in numerical models. Physically the roughness lengths 
correspond to the height above the surface where the assumed 
log profiles for the surface layer are valid. For momentum, 
there is a general dependence on wind speed such that the 
roughness length increases with increasing wind speed. Over 
water, the roughness elements consist of wave crests. This 
would seem to be easier to parameterize than inhomogeneous 
terrain over land, however as discussed in Fairall et al. 
(1996b) there are complicating factors such as changes in 
the wave structure in rough seas due to whitecapping and 
issues relating to wave age or swell direction when coming 
from a different direction than the surface winds. Assuming 
that Zo„=Zoh means that no distinction is made between the 

surface skin tenperature and the temperature at Field 

studies have shown that these values may differ by as much 
as 6 K over land. This temperature difference can be related 
to Zqj, and z^^ by the ratio which is given by Eq. (31) . As 

reported in Beljaars and Holtslag (1991) , for homogenous 

vegetated surfaces has been found to be about 6 which 
results in a z^JZ q^^-=\^ . This ratio is applied in the current 

version of COAMPS for over water grid points after 
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calculating the value of Zq^ using Charnock's relation and 


the u* from the previous time step. 




e. 



(31) 


Roughness lengths for momentum are usually derived from 
empirically derived climatological databases over land and 
can be calculated over water using Charnock's (1955) 
relation or other more recent formulations developed for 
light wind conditions as discussed in section II. The 
current TOGA-COARE parameterization uses an empirical set of 
equations that use the roughness Reynolds number (Rr) as 

input (Liu et al. 1979). These are based on laboratory 
measurements in wind tunnels, theoretical considerations 
about the interfacial siiblayer, and a requirement for 
relatively smooth transition between rough and smooth 
regimes. This relationship has been validated with a few 
field measurements as depicted in Figs. 1 and 2 but 
considerable scatter still exists. The application of a 
constant 1/10 ratio between momentum and heat roughness 
lengths results in the roughness lengths for heat continuing 
to increase in strong wind stress regimes (Rr>=0.8) in the 

Louis parameterization while decreasing in the TOGA-COARE 
method. 


26 




Figure 1. Normalized sensible heat roughness length (z 
denoted as Zoh in this thesis) versus'roughness Reynolds 
number (Rr) from model (line) and field measurements 
(circles). From Liu et al. (1979), 



Figure 2. Normalized latent heat roughness length (z^, 
denoted as Zq, in this thesis) versus roughness Reynolds 
number (Rr) from model (line) and field measurements 
(circles). From Liu et al. fl 975 l 
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III. METHODS 


A. THE MARINE AEROSOL CHARACTERIZATION EXPERIMENT (ACE-1) 

1. Overview 

The Southern Hemisphere Marine Aerosol Characterization 
Experiment (ACE-1) was conducted from 15 November to 14 
December 1995 in the South Pacific Ocean south of Tasmania. 
Coordinated measurements were made between the National 
Center for Atmospheric Research (NCAR) C-130 aircraft and a 
NOAA research vessel, R/V Discoverer. Details of the 
instrumentation and methodology are given in Wang et al. 
(1998a). Circular flight patterns were flown at four levels 
in the boundary layer and one level above. This was in order 
to provide longer duration at each altitude in order to 
improve confidence in the turbulence statistics and provide 
estimates of the large scale divergence from the mean wind. 
Some of the flights followed the southbound trajectory of 
tagged air columns so the flux observations crossed 
increasingly cold SST patterns and passed through regions of 
varying static stability and boiindary layer structure. Winds 
during the observation period were moderate to strong in the 
4-14 m s'^ range at a measurement height of 40 meters. 
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2 . 


Instmmentation 


Turbulence velocities were measured by the differential 
pressure technique using a gust probe mounted on the 
aircraft radome with a nominal accuracy of 0,15 m s'^ for 
vertical velocity and 1 m s‘^ for mean winds. A Rosemount 
platinum fast response temperature sensor was used to 
measure air temperature. Water vapor was measured with a 
Lyman-alpha hygrometer with a specified accuracy of ±4% for 
relative humidity and corrected using a chilled mirror 
hygrometer. SST measurements were made with a Heimann 
radiometer with a manufacturer specified accuracy of + 1 °C. 
Shortwave radiation was measured using an Eppley PSP 
pyranometer and longwave radiation using an Eppley PIR 
pyrgeometer. See NCAR (1995) for details of the aircraft 
instrumentation. 

B. DATA SELECTION 

The data analyzed for this study were chosen based on 
flights with a low level leg at approximately 40 meters and 
a relatively straight path or gradual turning along the leg. 
Flights in which the low level leg was at 100 meters or 
above were not used due to concerns that the measurements 
were possibly not made in the surface layer. In general the 
boundary layer height was observed to be between 500 and 
2000 meters during the ACE-1 period (Wang et al. 1998a). Our 
analysis based on Monin-Obukhov similarity theory required 
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measurements to remain within the surface layer, generally 
^ppi^oximated as the lowest 10% of the boundary layer. 

Flights with rapid course or speed changes were discarded 
because platform motion can produce error in the turbulence 
measurements. Table 1 is a summary of the mean conditions 
for the selected legs based on these criteria. The low level 
flight legs were then segmented into 28 km sections which 
overlapped by half of their length in order to examine the 
flux and mean variables across shorter distances. Bulk input 
variables were then calculated as the mean value for each 
shorter segment. This allowed us to examine mesoscale 
variations in the sea surface temperature, wind speed and 
static stability and their effects on the flux calculations. 
The aircraft speed was about 110 m s"^ so the segments 
correspond to about 5 minutes of flight time. 
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Table 1. Flight Data Summary. MSLP is Mean Sea-Level Pressure. 


Plight Leg 



Lon, 

Height 

(deg) 

(m) 



Wind 
Speed 
(m s*^) 



11 

01 

-53.5 

137.4 

40 

1013 

11 

05 

-53.2 

138.6 

38 

1012 

12 

01 

-48.2 

137.2 

39 

1027 

12 

06 

-48.5 

137.6 

41 

1027 

13 

01 

-49.4 

138.7 

44 

1000 

13 

06 

-48.4 

139.1 

61 

1001 

15 

01 

-47.7 

145.5 

39 

1014 

15 

06 

-47.1 

145.8 

40 

1015 

15 

07 

-46.8 

144.2 

39 

1016 

16 

01 

-54.1 

159.0 

39 

995 

16 

04 

-54.6 

158.9 

42 

995 

18 

04 

-45.0 

144.5 

39 

1004 

18 

00 

-45.2 

145.4 

40 

1006 

19 

01 

-46.2 

148.6 

37 

1003 

19 

06 

-46.8 

150.9 

45 

1002 

24 

01 

-45.2 

143.0 

41 

1010 

24 

06 

-45.2 

143.9 

41 

1011 

24 

09 

-45.2 

144.3 

42 

1011 

25 

01 

-45.9 

145.7 

38 

1010 


12.9 


12.9 


3.4 


4.5 


10.1 


13.2 



.3 


6.4 


10.4 


9.8 


.9 


10.3 


4.5 


4.8 



3.3 


3.8 


5.9 


6.2 


6.9 


6.5 


8.6 


8.9 


9.3 


4.2 


4.8 


11.4 


11.3 


11.1 


11.4 


10.8 


11.2 


11.4 


11.5 



2.0 


9.8 


9.6 


9.2 


9.2 


8.7 


10.1 


10.8 


10.1 
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C. THE EDDY CORRELATION METHOD 

Turbulent flvixes were calculated directly from aircraft 
data sampled at high frequency using the eddy correlation 
method. Statistical correlation is done on the perturbations 
or turbulent fluctuations of vertical velocity (w') , and 

wind (u' or v'), heat (6'), or moisture (g'). As reported in 

Wang et al (1998b), turbulence fluxes were calculated using 
25 Hz data measured by a C-130 aircraft. To calculate the 
turbulent fluxes, cospectra of the appropriate variables are 
obtained using the Fast Fourier Transform (FFT) technique. 
The cospectra are then integrated from the smallest 
resolvable scale (nominally 4 meters) to 6 km. The 6 km 
wavelength cutoff was selected to ensure all contributions 
from small-scale turbulence are included without 
incorporating mesoscale or larger scale variations. Figure 3 
is one example of the calculated cospectra. 

The length of the data segments used for spectral 
analysis was about 5 minutes which corresponds to 28 km sub¬ 
sections of the low level flight leg. For each circular leg, 
the adjacent 5 minute segments overlap by half of their 
length. In general, turbulent fluxes were found to vary 
considerably along each flight leg in the ACE-1 region as 
demonstrated in the example in Fig. 4. This variability 
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justifies the sectioning of the data into smaller regions in 
order to analyze the impact of the highly variable 
environmental conditions. Overlapping the segments and 
treating each subset as an independent datum produces a 
smoother transition for analysis and is appropriate for 
comparison to bulk values derived from the same process of 
overlapping sections. The corrected pressure altitude was 
about 40 meters for all flights and was considered to be in 
the surface layer. The boundary layer was complex and varied 
in height from 500 to 2000 meters in the region (Wang et al. 
1998b) but the assumption of surface layer constant flux is 
reasonable. 
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D. VARIATION OP SST AND THE WARM LAYER CORRECTION 

The Southern Ocean during the time and area of ACE-1 
contains several strong SST fronts and includes fairly 
complex oceanic mesoscale features. Due to the strong 
sensitivity of the bulk flux methods to the accuracy of the 
SST, it is essential to calibrate the SST to the best 
available accuracy. Initial comparison of the data with 
satellite-derived SST values indicated a possible systematic 
error in the aircraft data. One of the important tasks of 
this research study was therefore to obtain adequate 
corrections to the aircraft measured skin temperatures using 
all available measurements and analyses. 

In addition to the aircraft radiometric data, other 
sources of SST measurements in the ACE-1 region were 
satellite retrievals of Multi-Channel Sea Surface 
Temperatures (MCSST) from the NOAA-12 polar orbiter, 
conductivity/ temperature/ depth (CTD) casts by the R/V 
Southern Surveyor and the R/V Discoverer, and bulk inlet 
temperature measurements at a depth of 5 m from the R/V 
Discoverer. MCSST data provides the widest spatial and 
temporal coverage and is available for several of the 
flights used in this study. Figure 5 is a composite image of 
the satellite measured SST in the ACE-1 region. Since many 
of the features are masked by cloud cover at various times 
and the overall pattern changes very slowly, compositing the 
data collected during successive passes from 24 November to 
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8 December allows the major features to be seen. The strong 
spatial variation of SST is clearly depicted and the 
locations of the flights analyzed in this study are 
indicated by the text markers starting with "RF" and 
followed by the flight number. For the purpose of comparing 
satellite and aircraft SST measurements, precisely matching 
the location of the satellite data with the aircraft is not 
possible since the satellite data is an average temperature 
for a region several kilometers wide and the passes of the 
polar orbiter are several hours off from most of the low 
level flights. The CTD data has the best absolute accuracy 
and is closer to the surface than the inlet temperature. The 
CTD values may still not be exactly representative of the 
skin temperature due to the cool skin effect; however, the 
difference from the skin temperature is limited to a few 
tenths of a degree. The major disadvantage of the CTD data 
is that it was not collocated with the aircraft passes so 
exact comparisons were not possible. Inlet temperatures from 
the R/V Discoverer were recorded along the ship track every 
30 minutes. These values were measured with a thermistor 
with a manufactured specified accuracy of 0.001 ®C. The ship 
track overlapped the aircraft flight path on two occasions 
which makes direct comparisons between the platforms 
possible. 

The basic SST measurements used in this study were the 
C-130 radiometric skin temperatures since the aircraft data 
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provides variations along the flight track along which the 
surface flTOX parameterizations were examined in segments. 
Table 2 is a summary of the raw measurements from the 
various SST data sources. The satellite retrievals in Table 
2 were based on single satellite passes closest in time to 
the flight being analyzed and not on the 15 day average in 
Fig. 5. As indicated, the aircraft data consistently 
underestimates the surface temperature as compared to the 
satellite retrieved SST and ship inlet temperatures. This 
offset was corrected by closer analysis of the various data 
sources in order to approximate a constant correction factor 
to the C-130 data. The assumption here is that the offset in 
the C-130 data is constant for each flight. This appears to 
be reasonable since the variation in SST as depicted by the 
aircraft radiometer corresponds well with the variation in 
surface sensible heat flux calculated directly using the 
eddy correlation method. 

Due to the strong spatial variability of SST in the 
ACE-1 region, the ship inlet data should be ideal to 
calibrate the drift in the C-130 radiometric data in regions 
where the two platforms overlapped. This correction was done 
in Wang et al. (1998b) for the Lagrangian flights when the 

C-130 tracked one air column for two to three days. The R/V 
Discoverer and C-130 tracks were compared for the entire 
ACE-1 period and two days were identified as having 
sufficient overlap. As indicated in Table 2, the ship inlet 
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temperature was from 1,0 to 2.2 degrees higher than the 
aircraft radiometric temperature. Applying a +1.2 degree 
correction to the input SST for the bulk method indicated 
that the resulting sensible heat fluxes were still 
systematically underestimated when compared to the directly 
measured sensible heat flux. Increasing the SST correction 
further eliminated this discrepancy in the sensible heat 
flux, suggesting that some of the inlet temperature 
measurements underestimate the skin temperature value. This 
discrepancy was stronger when the inlet temperatures were 
collected during the late afternoon and in regions with weak 
wind stress. 

The finding that the ship inlet temperatures 
underestimated the skin temperature is expected. Price at 
al. (1986) showed considerable diurnal variation in the 
upper 20 meters of the ocean mixed layer. During the day, 
accumulation of the solar radiation in the upper mixed layer 
tends to form a stable warm layer in the upper few to twenty 
meters. If the inlet sensor is below the depth of the warm 
layer, the difference between the inlet temperature and the 
skin temperature can be significant and requires a 
correction. 

The warm layer correction from the TOGA-COARE algorithm 
provides for a means to adjust bulk inlet temperatures to 
more closely match required surface skin temperatures. This 
warm layer correction subroutine is based on a simplified 
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form of the Price-Weller-Pinkel turbulent mixing model 
(Price et al,, 1986). The warm layer correction is set to 
zero at local midnight and uses measurement time, downward 
solar and longwave heating, wind speed, sensor depth, and 
bulk water temperature to estimate the difference between 
the skin temperature and the temperature at measurement 
level. For wind speeds of 4-6 m s"\ the daytime surface was 
found to be on the order of 0.05-1.5 K warmer than the bulk 
temperature at 5 m and the layer depth was on the order of 
8-20 meters, with decreasing bulk-surface temperature 
difference and increasing layer depth values at increasing 
wind speeds (Fairall et al. 1996a). 

The magnitude of the warm layer effect was examined 
using the ACE-1 data. Solar and infrared radiative fliox 
inputs to the warm layer model were derived from aircraft 
measurements during flight 24 since the low level legs were 
in cloud free air and occurred during the main part of the 
local solar day. Wind stress, sensible heat and latent heat 
inputs were derived from aircraft data using the eddy 
correlation method for each flight where there was overlap 
between the ship and the aircraft. Similar corrections were 
also examined for a cloudy case by using observations from 
flight 12, which was in a region of moderate cloud cover. 

Figures 6(a)-6 (c) depict the temporal variation of the 
correction for three different environmental cases. Flight 
24 was in clear air with low wind stress, flight 15 was in 
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clear air with relatively high wind stress, and flight 12 
was in a cloud covered region with low wind stress. The 
turbulent fluxes used in the warm layer model were based on 
observed values for each flight. In addition to the diurnal 
evolution of the warm layer, the variation in the correction 
due to variation in the air-sea fluxes with the same diurnal 
heating cycle are considered and represented in each figure. 
The solid line represents the averaged values for the flight 
while the dashed lines were calculated using one positive or 
negative standard deviation of the observed momentum, latent 
heat, and sensible heat fluxes. The dashed lines therefore 
represent the likely upper and lower limits of the expected 
correction due to inhomogeneous turbulence flvix in the 
flight region. In general, the magnitude of the correction 
at this latitude and season is small, on the order of 0.2 K, 
except in the case of very light winds at a few hours past 
the maximum solar heating. Despite fairly opaque cloud cover 
during flight 12, flights 12 and 15 were found to have 
fairly similar corrections under similar wind stress 
conditions. Insolation was reduced by about 40% during 
flight 12, however downwelling longwave radiation was 
stronger due to the emittance of the cloud layer resulting 
in a net heating effect on the ocean mixed layer similar to 
the clear air case. The warm layer effect was much larger 
during flight 24 with the same clear conditions. This was 
due to weaker turbulence mixing due to wind stress resulting 
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in a shallower layer. Because of the strong non-linear 
response of the warm layer model to wind stress and the high 
variability of wind stress along the flight legs, a 
hypothetical case was run using averaged latent and sensible 
heat fluxes and varying the wind stress by 0.01 N m'^ 
increments. As can be seen in Fig. 6(d), the warm layer 
effect increased dramatically for wind stresses below 0.02 N 
m , which is equivalent to wind speeds on the order of 4 m 
s"^. 

The maximum warm layer correction depicted in Fig. 6(d) 
varied from 0.3-1.3 K due to different magnitudes of the 
specified wind stress. The warm layer corrections listed in 
Table 2 for flights 15 and 24 are based on the local solar 
time of the shipboard observations and the mean wind stress 
for the segments of the flight leg near the ship location. 
These values should be used with caution due to uncertainty 
arising from the inhomogeneity of the turbulent fluxes and 
the temporal evolution at the shipboard measurement 
location. 

Based on several independent sources of SST data 
discussed above and the observed linear offset in the 
sensible heat flux values between the bulk methods and the 
eddy correlation data, a correction of +3.2 °C was applied 
to the C-130 measured SST for flight 11 and a correction of 
+2.2 °C was applied to flights 12-25. 
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Sea surface temperature composite 24 Nov - 8 Dec 95 

ACE-1 cruise tracks. Discoverer: white line, Southern Surveyor: black line, aircraftidotted lines 

__||^ ^ 145 



Figure 5. Composite MCSST imagery of the mesoscale oceanic variability 
during ACE-1. Adapted from Griffiths et. al (1998). 


44 
















E. THE BULK FLUX CALCULATIONS 

The high-rate data collected from the C-130 was reduced 
to 1 Hz data by averaging the values for each second. The 
low rate data for the overlapping 28 km segments of the low 
level flight legs was then averaged and these values were 
used as input to the bulk methods. The calculated fluxes are 
output at measurement level which is 40-60 meters. Transfer 
functions are brought down to a standard level of 10 meters 
and neutral conditions for comparison purposes. 


Table 2. Sea Surface Temperature (®C) from various platforms and 
calculations. 


Flight 

Lat 

Lon 

C-130 

Satellite 

Ship 

Warm Layer 
Corr. 

11 

-53 

138 

1.1 

4.5 

— 

— « * 

12 

-48 

137 

6.4 

8.5 

- - ^ 

... 

13 

-49 

139 

6.6 

9.0 

— 

... 

15 

-47 

146 

6.8 

,9.5 

8.8 

0.1 

16 

-54 

159 

1.7 

4.0 

— 

... 

18 

-45 

145 

9.7 

12.0 

— 

> ... 

19 

-46 

150 

9.2 

12.0 

. . ^ 

... 

24 

-45 

143 

8.7 

11.0 

10.9 

0.4 

24 

-45 

144 

10.4 

12.0 

11.4 

1.1 

25 

-46 

146 

10.1 

12.5 

— 

— 
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IV. 


RESULTS 


The data collected during the ACE-1 Intensive 
Observation Period (lOP) covers a wide range of conditions 
for static stability and wind speed. The highest fluxes 
occurred during flights 11, 12, and 13. Flights 11 and 13 
are closer to neutral static stability with moderately high 
wind speeds while flight 12 is strongly convective with 
relatively light winds. Table 3 is a summary of mean 
conditions and air-sea fluxes for the low level legs 
examined. The listed fluxes are derived from the eddy 
correlation method and are not corrected to a standard level 
but should be close to the values at 10 meters since they 
were measured in the surface layer. The data in Table 3 is 
an average value for the entire 180 km leg while the data 
used in the scatter plots in Figs, 10-15 discussed below is 
segmented into overlapping 28 km lengths as described in 
section III. 

A. TURBULENT FLUXES 

Figures 7-9 depict the momentum, sensible heat, and 
latent heat fluxes for all flights examined. The Louis and 
COARE methods follow the trend of the variability in the 
fluxes in this region quite well but both show error in the 
absolute magnitude. As seen in Fig. 7, flights 11 and 13 
were in regions of high wind stress. The bulk methods 
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produced the largest discrepancies of the wind stress in 
these flights. The COARE algorithm comes closer but fails to 
capture the entire range of values. Figure 8 depicts a 
similar trend for sensible heat flux as that seen for wind 
stress. Neither method captures the maximum values seen in 
flight 11 and 13 however the COARE method comes closer. 
Figure 9 is a comparison of the latent heat flux. Flight 13 
again shows large fluxes that are underestimated by the bulk 
methods but the Louis and COARE estimates are much closer to 
each other than to the observed values. Moderate latent heat 
fluxes are common over more flights than the sensible heat 
and momentum fluxes. Latent heat fluxes are overestimated in 
flights 18, 19 and 24 by both bulk methods however the Louis 
method comes closer to the observed values. 

The scatter plots in Figs. 10-12 are direct comparisons 
of the calculated turbulent fluxes against the observed 
values and between the two bulk estimates. As seen in Fig. 

10(a), in high wind stress conditions the Louis method 
underestimates the stress by about 50%. The COARE method 
depicted in Fig. 10(b) does slightly better but still shows 
large errors. Comparison of the COARE and Louis estimates in 
Fig. 10(c) shows that they agree very closely in low stress 
regions and begin to differ when the momentum flux is above 
0.2 N m”^. This corresponds to flights 11 and 13 where mean 
winds were greater than 11 m s'^. 

As seen in Fig. 11, in high sensible heat fl\ix 
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conditions the bulk methods underestimate the heat flxix by 
about 30%. The COARE method does slightly better (Fig. 
11(b)), especially with fliixes above 15 W m"^, and produces 
larger magnitude negative fluxes in the stable regime. In 
low heat flux conditions, there is considerable scatter in 
the data, with COARE doing better in some cases such as 
flight 19, and Louis doing better in others such as flight 
15 and 16. There are not enough cases of stable static 
stability in this dataset to produce a definitive answer but 
the COARE estimates do seem to verify better. 

In Fig. 12, the Louis method clearly does better in 
estimating the latent heat flux. Again neither method 
captures the high fluxes in flight 13 and latent heat flux. 
estimates show more scatter than wind stress or sensible 
heat flux. For fluxes less than 75 W m'^, the COARE method 
overestimates the latent heat flux by about 30-40%. 

Figures 13-15 are comparisons of the calculated fluxes 
as a function of wind speed and static stability. As can be 
seen in Fig. 13(a), the regions of high wind stress were 
neutral to weakly convective. The results in Figs. 13-15 
indicate that the turbulent flux dependence on stability can 
be separated into two regimes: the high-wind regime and the 
low-wind regime, although in the high-wind regime the 
stability is limited to close to neutral. These results 
appear to indicate that the resultant stress is a strong 
function of z/L in the high-wind conditions, although more 


49 



observations are needed for a definitive conclusion. In the 
low-wind regime, however, the stress appears to be 
insensitive to variations in static stability. In Figs. 14 
and 15, large fluxes of heat and moisture are seen in near 
neutral conditions with reduced values in moderately 
convective conditions and increasing values in more 
convective regimes. These heat fluxes can also be separated 
into two groups based on low and high wind conditions with 
the division between these two regimes nominally at 10 m s' 
The high-wind regime is where the calculated and observed 
fluxes from the three methods differs the most. The wind 
stress as computed by the bulk formulae is roughly 
proportional to U^, which would indicate nearly constant 
As seen in Fig. 13(b), this behavior is not seen in the eddy 
stresses which increase much faster with wind speed. 

Sensible heat flux appears to decrease with wind speed for 
speeds less than 10 m s'^ and then increase with wind speed 
above this value. Latent heat fliox seems to show a similar 
pattern to the sensible heat flux. While considerable 
scatter was seen in the actual wind stress, the bulk methods 
correlated very closely and fit the trend of the eddy 
correlation data well. In Fig. 14, the sensible heat flux 
increased rapidly above 11 m s'^, despite near neutral 
stability, which is the opposite of what was found in the 
HEXOS data (De Cosmo et al. 1996). Figure 15 shows a modest 
increase in latent heat flux in stronger wind conditions as 
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well as with strongly convective conditions in light winds. 


B. TRANSFER COEFFICIENTS 

The bulk transfer coefficients at a standard height of 
10 meters are depicted in Fig. 16. Due to the day-to-day 
variation of the environmental conditions such as static 
stability and wind speed, we observe large variation in the 
transfer coefficients from flight to flight. However, the 
COARE transfer coefficients do show more variability from 
flight to flight and within each flight than the Louis 
values. However, this does not necessarily suggest that the 
COARE method respond faster to environmental changes than 
the Louis method as both methods depict similar variations 
in flux as seen in Figs. 7-10. As expected from the scatter 
plots of resultant fluxes in Figs. 11 and 12, the COARE 
method produces higher latent and sensible heat transfer 
coefficients than the Louis method for most of the observed 
cases but over-estimates the latent heat transfer 
coefficient and under-estimates the sensible heat 
coefficients in some regimes. The bulk methods also do not 
capture the variability represented in the direct flux 
measurements. As can be seen in Figs. 17 and 18, the flux 
transfer coefficients are found to be a weak function of 
stability except in moderate to strong wind conditions, 
although at these wind speeds the stability is restricted to 
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near neutral. However, the latent and sensible heat transfer 
coefficients are seen to increase with increasing wind 
speed. The general trend of two distinct regimes for 
sensible and latent heat flux is represented by the eddy 
correlation method as well as the bulk methods. The high 
wind, near neutral case and the low wind, highly convective 
case both result in required flux transfer coefficients 
based on the eddy correlation method that are much higher 
than those produced by either bulk method. In between these 
two regimes the formulation by Louis matches the 
observations better for latent heat flux (Fig. 17(c)) while 
the COARE formulation matches the sensible heat flux (Fig. 
17(b)) more closely. As can be seen in Fig. 17(b) and 17(c), 
the thermal transfer coefficients in the bulk methods do not 
increase rapidly enough in the highly convective regime for 
z/L values less than -0.5. 

C. ROUGHNESS LENGTHS 

Figure 19 compares the roughness lengths for momentum, 
heat and moisture for the Louis and COARE methods. As seen 
from Fig. 19(a), the Louis method calculates slightly 
smaller values for momentum roughness length than the COARE 
method. This is expected since the COARE algorithm adds an 
additional term from the LKB formulation for smooth flow 
with the difference being the greatest at low wind speeds. 
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For heat (Fig. 19(b)) and moisture (Fig. 19(c)) roughness 
lengths, the differences are much larger as expected since 
the Louis method uses a fixed fraction of momentum roughness 
length to calculate latent and sensible heat roughness 
lengths while the COARE algorithm uses an empirically based 
parameter!zation. 

D. DISCUSSION 

It is noted in Figs. 18(b) and 18(c) that the 
variations of and Cg with wind speed are similar to the 
variations of Zoh and Zoq with Reynolds roughness number (Rr) 
in Figs. 1 and 2. This suggests that the observed 
differences in the transfer coefficients and resultant 
fluxes may result from the different formulations of the 
roughness length in the two methods. The Louis 
parameterization was thus modified to calculate the 
roughness length using the COARE formulations in Eq. (21) 
and Figs. 1 and 2. Figure 20 is a comparison of the COARE 
roughness lengths with the roughness lengths calculated from 
this modified version of the Louis method. The modified 
Louis values are much closer to the COARE values and depict 
more variability than the original Louis values. They still 
do not match the COARE roughness lengths exactly because the 
initial value for u» is different between the two methods 

but the momentum roughness length in Fig. 20(a) is much 
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closer at the low values where the second term from the LKB 
method has the most effect. 

Figures 21-23 compare the resultant fluxes from the 
modified Louis to the COARE and eddy correlation data. The 
modification of momentum roughness length has little impact 
on wind stress estimates between the two methods for the 
conditions observed the ACE-1 region. As can be seen from 
Fig. 21(c), the wind stress from the two bulk methods 
matches slightly more closely at low values although the 
Louis method still has lower values above 0.2 N m‘^. For 
sensible heat flux (Fig. 22) the Louis method still under 
estimates the flux at high values. The modified sensible 
heat roughness lengths do cause the Louis estimates to match 
the COARE results more closely. Modifying the thermal 
roughness length does not seem to improve the observed 
underestimate of the Louis method for the stable conditions. 
For latent heat flux, the parameterization in the COARE 
algorithm for Zoq may have a significant negative impact on 

the COARE results. This is supported by comparing Fig. 12 
with the results in Fig. 23. The Louis latent heat fluxes 
are very similar to the COARE results when the same 
formulation for Zoq is used. We therefore conclude that the 
fact that the COARE algorithm overestimates the latent heat 
flux is likely caused by am inappropriate moisture roughness 
length (Zoq) parameterization. 
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Table 3. Summary of average aircraft loaction, stability, wind 
conditions, and Air-Sea Fluxes. z/L is an estimate of stability 
based on measurement height (z) and Monin-Obukhov Length (L) as 
defined in Eg. (7) ._ 


Fit. 

Lat. 

(deg) 

Lon. 

(deg) 

Wind 
Speed 
(m s'^) 

z/L 

IqH 

Tau . 

(N m-*) 

H 

(W m"^) 

E 

(W m"2) 

11 

-53.4 

137.8 

12.9 

-0.05 

-0.71 

0.27 

13.78 

56.23 

12 

-48.3 

137.4 

4.0 

-3.11 

-2.14 

0.06 

16.78 

53.18 

13 

-48.9 

138.9 

11.6 

-0.20 

-1.40 

0.30 

32.56 

105.59 

15 

-47.4 

145.6 

5.5 

0.01 

-0.18 

0.05 

2.33 

26.82 

16 

-54.1 

159.0 

7.2 

0.12 

0.08 

0.04 

8.22 

37.92 

18 

-45.1 

145.0 

10.1 

-0.05 

-0.40 

0.13 

2.07 

48.09 

19 

-46.5 

149.8 

9.2 

0.04 

-0.04 

0.10 

-1.53 

28.27 

24 

-45.2 

143.7 

4.6 

-0.82 

-0.88 

0.03 

2.69 

21.28 

25 

-45.9 

145.7 

4.9 

-0.58 

-0.73 

0.03 

2.16 

22.37 
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Figure 7. Calculated and observed wind stress vs. flight 
segment for all analyzed flights during the ACE-1 lOP. 
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Figure 10. Comparison of the calculated wind stress (N m'^) 
using two bulk methods and the observed flux from the eddy- 
correlation method. 
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Figure 11. Same as Fig. 10, except for sensible heat flux. 
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Wind Stress vs Wind Speed 



10 meter Wind Speed 


Figure 13. Variations of wind stress as a function of (a), 
static stability and (b), 10 m wind speed, for all analyzed 
flights during the ACE-1 lOP. The method used to obtain the 
fluxes is indicated in the legend. 
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Figure 14. Same as Fig. 13 except for sensible heat flux. 
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Sensible Heat Drag Coefficient 



Figure 16. Bulk transfer coefficients at 10m height obtained 
from the Louis and COARE methods analyzed from flights 
during the ACE-1 lOP. Flight numbers are indicated in (a). 



















10m Stability z/L 


Figure 17. Bulk transfer coefficients at 10m height as a 
function of static stability. The method used to obtain the 
transfer coefficients is indicated in the legend. 
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Figure 18. Neutral bulk transfer coefficients at 10m height 
as a function of 10m wind speed. The method used to obtain 
the transfer coefficients is indicated in the legend. 
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Figure 19. Comparison of roughness lengths for momentum 
(zoro) , heat (zoh) and moisture (zoq) between the Louis and 
COARE algorithms. 
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Figure 20. Comparison of roughness lengths for momentum 
(zo) , heat (zqh) and moisture (zoq) between the Louis method 
modified to use the COARE parameterization and from the 
COARE algorithm. Differences are due to different initial 
estimates of u. in the two methods. 
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Figure 21. Comparison between the calculated wind stress (N 
m'^) using the modified Louis method and the COARE algorithm 
and the observed wind stress from the eddy correlation 
method. 
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Figure 22. Same as Fig. 21, except for sensible heat flxuc. 
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Modified Louis Latent Heat 




Figure 23. Same as Fig. 21, except for latent heat flux. 
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V. 


SUMMARY AND CONCLUSIONS 


The Louis et al. (1982) fliix parameterization as 
applied in the COAMPS model was compared to the TOGA-COARE 
flux algorithm using observations collected during ACE-1, 15 
November 1995 to 15 December 1995 in the South Pacific Ocean 
southwest of Tasmania. The Louis method is computationally 
very compact and works well for previous forecast models. 
However, improvements in numerical models and the increasing 
effort in modelling coupled air-sea systems requires 
improved accuracy in the air-sea flux calculations. The 
TOGA-COARE flux algorithm is a state-of-the-art method 
incorporating several physical considerations neglected in 
earlier turbulent fl\ix parameterizations. This method is 
specifically designed to produce more accurate flux 
estimates from bulk methods but is an iterative approach 
that is more computationally expensive. As this method is to 
be incorporated into the U. S. Navy's Coupled Ocean- 
Atmosphere Mesoscale Prediction System (COAMPS), it is 
desirable to evaluate its performance relative to the 
existing fl\ix parameterization (Louis at al., 1982) and 
direct observations. 

These two methods were found to be extremely sensitive 
to accurate sea surface temperature (SST) measurements. A 
detailed analysis of all available SST data was conducted in 
order to ensure adequate accuracy of the specified skin 
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temperature from the observations. Comparison of shipboard 
and satellite-derived SST estimates with the aircraft data 
indicated that the aircraft data had a consistent low bias. 
The warm-layer correction algorithms incorporated in the 
TOGA-COARE method were used to adjust the observed ship 
inlet temperatures to skin temperature for this comparison. 
The modeled warm layer correction based on observed 
radiative warming and surface turbulent fluxes was found to 
vary from 0.1 to 1.3 °C in the ACE-1 region. This 
correction was particularly sensitive to wind stress and was 
minimal for wind stresses above 0.03 N m"^. The resultant 
adjustment of the aircraft observations of sea surface 
temperatures was found to greatly improve the resultant 
fluxes. 

The two bulk methods showed close agreement with each 
other and with eddy correlation measurements of the actual 
fluxes for momentum flux. For latent and sensible heat flux, 
the two methods showed good agreement in low flux regimes 
but both methods underestimated the flux for values above 20 
W m"^ for sensible heat flux and 70 W m'^ for latent heat 
flux. 

The high heat and moisture flux cases can be separated 
into two regimes, cases with wind speeds above 10 m s"^, 
which were also characterized by near neutral stability, and 
cases with low wind speeds, which were also highly 
convective. Calculation of the required transfer 
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coefficients from the directly observed fluxes indicates 
that the transfer coefficients in the bulk methods do not 
increase rapidly enough for z/L values less than -0.5. In 
the high wind regime, the transfer coefficients are seen to 
be very sensitive to static stability. 

In general, the COARE method was found to be better 
than the Louis method in estimating the sensible heat fl;ax 
while the Louis method was superior for latent heat flux. 

The observed decrease in latent and sensible heat transfer 
coefficients with wind speed in the COARE model was similar 
to the decrease in heat and moisture roughness lengths with 
u* in the COARE parameterization. The Louis method was 

modified to match the more sophisticated COARE 
parameterization for momentum, thermal, and moisture 
roughness lengths. Little improvement was found for 
estimates of wind stress, which is expected since the 
additional term in the COARE model accounts for smooth flow 
regimes which were not observed in this study. Sensible heat 
flux estimates were observed to improve slightly in the 
modified Louis method, especially at higher flxix levels, but 
modifying the thermal roughness length in the Louis method 
did not improve the sensible heat flux estimates in stable 
conditions. The modified Louis method was found to 
overestimate the latent heat flux in the same manner as the 
COARE algorithm which indicates that the parameterization 
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for latent heat roughness length in the COARE algorithm 
requires further study. 

The observations and analysis in this study indicate 
that bulk flux parameterizations seem to break down in high 
wind and highly convective conditions. This is not 
surprising as the field studies and theoretical treatments 
used in developing Monin-Obukhov similarity theory are often 
based on mean synoptic environments and surface layer 
profiles. Accurate estimates of heat and moisture flux in 
these regimes are particularly sensitive to the correct 
formulation for heat and moisture roughness lengths and our 
observations indicate that moisture roughness length in 
particular requires further study. In high wind conditions, 
static stability approaches near neutral conditions, however 
heat and moisture fluxes are still quite strong, and an 
accurate estimate of these fluxes is very sensitive to small 
variations in the stability parameter (z/L). This in turn 
requires very accurate specification of the air-sea 
temperature contrast, which requires better SST measurements 
than are currently available in operational models. 
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